Mon. Not. R. Astron. Soc. 000, EM^J (2003) 



Printed 2 February 2008 



(MN BT E X style file v2.2) 



Outer-magnetospheric model for 
Formation of sub-GeV spectrum 



Vela- like pulsars: 



J. Takata, 1 * S. Shibata, 2 K. Hirotani, 3 



1 Graduate School of Science and Engineering, Yamagata University, Yamagata 990-8560, Japan 



2 Department of Physics, Yamagata University, Yamagata 990-8560, Japan 



3 Max- Planck- Institut fuer Kernphysik, Postfach 103980 D-69029 Heidelberg, Germany 



Released 2003 Xxxxx XX 



ABSTRACT 



We investigate the 7-ray emission from a outer-gap accelerator, which is lo- 
cated in the outer-magnetosphere of a pulsar. The charge depletion from the 
Goldreich- Julian density causes a large electric field along the magnetic field 
lines. The electric field accelerates electrons and positrons to high energies. 
We solve the formation of the electric field self-consistently with the radia- 
tion and the pair-creation processes in one-dimensional geometry along the 
magnetic field lines, and calculate curvature spectrum to compare with the 
observations. 

We find that because the particles escape from the accelerating region 
with large Lorentz factors as 10 7 5 , curvature radiation from the outside of 
the gap, in which there is no electric field, still has a important contribu- 
tion to the spectrum. Including the outside emission, the model spectra are 
improved between 100 MeV and 1 GeV for the Vela-type pulsars. Since the 
field line curvature close to the light cylinder affects the spectrum in several 
hundred MeV, there is a possibility to diagnose the field structure by using 
the spectrum in these bands. The difference found in the 7-ray spectra of 
the twin pulsars, PSR B0833-45 (Vela) and B1706-44, is deduced from the 
dimensionless electric current running through the gap. The extension to the 
three-dimensional modeling with self-consistent electrodynamics is necessary 
to understand the difference in the observed flux between the twin pulsars. 



Key words: gamma-rays: theory -pulsars: individual (Vela pulsar, PSR B1706- 
44). 
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1 INTRODUCTION 

The Compton Gamma Ray Observatory had detected seven 7-ray pulsars (Thompson et al. 
1999). The observed light curves and energy spectra tell about the energy of accelerated 
particles and the geometry of 7-ray emission regions in the pulsar magnetosphere. These 
observations provide constraints on proposed radiation models. The next-generation 7-ray 
space telescopes and ground based Cerenkov telescopes are expected to provide further 
constraints on the models. 

A pulsation in 7-rays implies that the acceleration of particles and the subsequent radi- 
ation occur within the light cylinder, of which axial distance is given by zui c = c/Q, where 
Q is angular frequency of the star, and c the speed of light. 

The charge depletion from Goldreich- Julian charge density causes a large electric field 
Eu along the magnetic field lines. The Goldreich- Julian density is given by (Goldreich & 
Julian 1969, here after GJ) 
QB Z 1 

PGJ 2txc [1 - (ttw/c) 2 ] ' {> 
where B z is the component of the magnetic field along the rotation axis, and w the axial 
distance. A potential drop in such a charge depletion region (gap) is only a small fraction 
of the available electromotive force exerted on the spinning neutron star surface. The 7-ray 
emissions in the pulsar magnetosphere have been argued with polar-cap model (e.g. Sturrock 
1971; Ruderman & Sutherland 1975) and outer-gap model (e.g. Cheng, Ho & Ruderman 
1986a, b); each acceleration region is located near stellar surface above the magnetic poles 
and in the outer-magnetosphere around the null charge surface (B z = 0) above the last-open 
field line, respectively. 

After Cheng et al. (1986a,b), the outer-gap model has been studied in two different 
pictures. In the first picture, the observed light curves have been studied in an effort to con- 
struct the thee-dimensional models for the outer-gap (Romani & Yadigaroglu 1995; Romani 
1996; Cheng, Ruderman & Zhang 2000). The peaks in the light curves are interpreted as an 
effect of aberration and time delay for emitted photons, thereby the outer-gap model has 
been successful in explaining the double peaks in one period and phase separation of the 
pulses in different energies from the radio to 7-ray. The observed phase resolved spectra have 
also been calculated with the models, in which pair cascade process has been solved with 
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assumed field-aligned electric field E\\. The predicted phase resolved spectrum depends on 
the assumption of the strength and distribution of En. Thus, the models are not satisfactory 
in the sense of electrodynamics. 

The other approach has focused on electrodynamics in the outer-gap. Cheng et si. 
(1986a) discussed the electric structure in the vacuum gap. They gave a solution for the 
Poisson equation with rectilinear geometry. Hirotani & Shibata (1999, hereafter HS) pro- 
posed a non-vacuum electrodynamical model for the outer-gap. Although they worked on 
only one dimension along the last-open field line, they solved the accelerating electric field 
self-consistently with curvature radiation and pair-creation processes. Although the calcu- 
lated 7-ray spectra were consistent with the EGRET observation in the GeV band, there 
were much harder than what is observed below 1 GeV. 

In this paper, to improve the HS model, we take account of the curvature radiation from 
the outside of the gap. Even though the electric field is screened out in this region, particles 
still possess sufficient energy and hence contribute to the 7-ray luminosity between 100 MeV 
and 1 GeV. Quite recently, Hirotani, Harding & Shibata (2003) presented a model in which 
particle motion is precisely solved under the electric force and radiation reaction force (see 
also Hirotani 2003 for a review). Although they also took account of the radiation from the 
outside of the gap, we treat this effect in some detail by considering non-dipole field near the 
light cylinder because the field line are expected to deviate from dipole configuration by the 
influence of the pulsar wind or the electric current. Then, we show that the position of low 
energy cut-off in the curvature spectrum depends on the curvature radius of the field lines 
near the light cylinder. In applying our model, we consider the twin pulsars, PSR B0833-45 
(Vela) and B1706-44. It is remarkable that these pulsars show distinctive 7-ray spectra even 
though they have very similar magnetospheric parameters as well as surface temperatures 
(Gotthelf, Halpern & Dodson 2002). Thus this pair of pulsars gives unique opportunity to 
examine the model capability. We shall show that the difference of spectral peak energy 
between the two, about 3 GeV for Vela and about 1 GeV for B1706-44, is understood well 
with the present model. On the other hand, we elucidate limitations of one-dimensional 
approximation. 

In £0we describe the basic equations for the electrodynamics in the gap. In £0we discuss 
the curvature radiation from the outside of the gap. In §3] we apply our model to the Vela 
pulsar and PSR B1706-44. The necessity for the three-dimensional electrodynamical model 
is discussed in £0 
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2 ONE DIMENSIONAL OUTER-GAP MODEL 



In this section, we present the basic equations to describe the electrodynamics in the gap, 
after HS. We give the Poisson equation for the acceleration field (§ §2.1JI . the continuity 
equations for particles and 7-rays (§ §2.2J) . We impose the boundary conditions in § §2.31 and 
finally describe the X-ray field in § §2.41 



2.1 Poisson equation 

If the magnetosphere is corotating with the star, an electric field, 
(n x r) x B „ x 

E = -± j = -V$ co , (2) 

which is perpendicular to the magnetic field B, is exerted at position r. The corotational 
potential <3? co satisfies the Poisson equation 

V 2 $ co = -4vrp GJ . (3) 

Meanwhile in the gap, since the real charge density p deviates from pcj, the electrostatic 
potential $ differs from $ co . The difference is extracted as the non- corotational potential, 
&nco = $ — $co- The Poisson equation, V 2 $ = — 47rp, yields 

V 2 $ nco = -47rp e// , (4) 

where p e ff = p — pcj is the effective charge density which produces the field-aligned electric 
field, £,| = —B ■ V$ncol\B\. 

If the gap width W\\ along the magnetic field lines is much shorter than the curvature 
radius (~ wi c ) of the field lines, we may assume that the field lines are straight in the gap. 
In §H we shall see that W\\ is about 10 per cent of w\ c . By assuming that the gap width is 
much smaller than the trans- field thickness of the gap, we rewrite equation (JIJ) as 

n.nn \ o 1 



-4tt [p(s) - p G j(s)} 



ds 2 

PGj{s) 



-Aire 



N + (s)-N(s)- 



(5) 



where N + (iV_) is the positrons (electrons) number density, s is the arc length along the 
last-open field line from the surface, and — e is the electron charge. The accelerating electric 
field is given by 
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Unless the gap is close to the light cylinder, the GJ charge density (JTJ) is approximated 

as 

Pgj - -2^. (7) 
We use this expression for the GJ charge density. 

2.2 Continuity equations for particles and 7-rays 

As the Poisson equation © indicates, particle's distribution is needed to obtain the electric 
structure in the gap. If charged particles move along the magnetic field lines with the speed 
~ c, the steady continuity equation without any particle sources yields 

constant along magnetic field line. (8) 



N ± (s 



B(s) 

However, 7-rays radiated by particles may convert into electron-positron pairs. In the outer- 
magnetosphere, X7 — > e~e + process has a particular importance with X-ray illumination 
from the neutron star surface or the magnetosphere. Taking account of the pair-creation, we 
write down the continuity equation into the form, 

±B ^(§) = ^*f^^ G++ *- G -i' < 9) 

where we assume that the electric field is positive in the gap for the definite of sign of 
polarity so that e + (or e - ) may be accelerated outward (or inward); G p+ (s, e 7 ) (or G P J) 
denotes the distribution function of outwardly (or inwardly) propagating 7-rays; e 7 is the 
photon energy normalized by m e c 2 ; ^> is the angle between the particle's motion and the 
meridional plane and is assumed to be \I/ = sin _1 (rf2sin#/c), where 6 is the colatitude at 
the point considered. The pair-creation rate %,±(s,e 7 ) per unit time per 7-ray photon with 



e 7 is written down as 



f°° dl\x 

V P ±(s,e 1 ) = (1 -/i c )c / de x —, — a p , (10) 

Jtth <lex 

where dex • dNx/dex is the X-ray number density between energies m e c 2 ex and m e c 2 (ex + 
dex) , cos -1 fi c is the collision angle between a X-ray and a 7-ray, m e c 2 e t h is the threshold 
energy for the pair-creation, and o~ p is the pair-creation cross-section, which is given by 

1-1-9) _ 1 

(11) 



3 

o- p (e y ,e x ) = —Ml - 1,2 



(3-77 4 )ln^J -2v(2-v 2 ) 



u(e 7 ,e x ) = Jl - — — , 

V 1 - fJ, c e 7 ex 

where o~t is the Thomson cross section. 
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Let us next consider the continuity equation for 7-rays. The 7-rays radiated by high en- 
ergy particles are beamed nearly tangentially to the magnetic field. Although 7-rays deviate 
from the magnetic field line on which 7-rays were emitted, the one-dimensional approxima- 
tion does not treat this effect correctly. We deal with 7-ray distribution along the magnetic 
field line. Then we can write down the continuity equation into the form, 



±B{s)±- 
as 



G±(s,e 7 ) 
B(s) 



r] p ±G±(s, e 7 ) + r] c N±(s) ^ 



ccos ^> 

where the first term of the right-hand side represents annihilation of 7-rays by the pair- 
creation, and the second term represents production by the curvature radiation. The emis- 
sivity for the curvature radiation r] c is given by 

Vc(R c ,T,e,) = ^^F(^), (13) 



hR c e 7 



3 



3 hcF 

47r m e c 2 R c 



(14) 



and 



poo 

F{x)=x K 5/3 (t)dt, (15) 

J X 

where R c is the curvature radius of the magnetic field line, T is the Lorentz factor of particles, 
K 5 / 3 is the modified Bessel function of order 5/3, h is the Planck constant, and m e c 2 e c gives 
the characteristic curvature photon energy radiated by a particle with Lorentz factor T. The 
Lorentz factor in the gap is obtained by assuming that the particle's motion immediately 
saturates in the balance between the electric and the radiation reaction forces, 

/od2 \V 4 

T sat (R c ,E ll )=i-^E ll + l\ . (16) 

The unity in the brackets in (|16|) does not influence T sat in the gap except for near the 
boundaries, where En =0. 

Saturation simplifies the problem significantly, but requires a condition. The r sot is the 
Lorentz factor at which the typical accelerating time t ac equals the typical traditional damp- 
ing time td of the particles. The particles whose Lorentz factor is lager than T sat will be 
decelerated by the radiation back reaction because td < t ac . On the other hand, the particles 
which have T < T sat in the gap will be accelerated because td > t ac . If the two timescales 
are less than crossing time t cr = W\\/c, 

tar > t ac and t cr > t d , (17) 
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then the particles in the gap will have the saturated Lorentz factor given by equation (|16|). As 
has been recently showed by Hirotani et al. (2003), in which motion of unsaturated particles 
are solved, the condition t cr 3> t ac in (fTTj) is not satisfied effectively for some pulsars and/or 
part of the gap. In the present application to the Vela-type pulsars, 7-rays producing the 
peak in the spectrum emerge from the central part of the gap, for which use of (j!6|) provides 
a good approximation. However, emission near and the outside of the gap requires a different 
treatment, which is discussed in JJJ 

2.3 Boundary conditions 

In this subsection, we give the boundary conditions to solve equations ©, ©, and (fT2"j) . 
The inner (s — si) and outer (s = S2) boundaries are defined so that E\\ may vanish, 

S||( Sl ) = S||(s 2 ) = 0. (18) 

We assume that 7-rays do not come into the gap across either of the boundaries, that is, 

G + ( Sl ) = G_( S2 ) = 0. (19) 

Attempting to explain the activity of pulsars, it is generally believed that some currents 
circulate in the magnetosphere globally. Hence we permit particles to come into gap through 
the boundaries. We parameterize the inflows across the inner and the outer boundaries such 
that, 

»+M ju »^ =j ,. (20 ) 



flB(si)/2vrc J ' OB(s 2 )/27rc 
If ji = 1 (or j 2 = 1), particles are injected into the gap at the Goldreich- Julian rate at s — S\ 

(or s = S2). The particle's continuity equations Q give the current conservation per unit 
flux tube, 

N+(s) NJs) 

hot = — r~r-, 1 - — r~r-, = constant along s. (21) 

By defining j g as the amount of created current carriers per unit flux tube within the gap, 
the conservation low (|2*T|) yields 

jl + 32 + jg = jtot- (22) 

The global conditions, which include an interaction with the pulsar wind, should determine 
the current. In our local model, therefore, we adopt (jtot, ji, 32) or (j g , 31,32) as the set of 
free model parameters. 

We have the 7 boundary conditions (|18j) ~ (|22|) against the 5 unknown functions. This is 
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because E\\ = is imposed at both boundaries and because j tot (or j g ) is externally imposed. 
In other words, we cannot impose the conditions E\\(si) = and E\\(s2) = for arbitrary 
given boundaries si and S2- By moving the boundaries step by step iteratively, we seek for 
the boundaries that satisfy the required conditions. 

2.4 X-ray field in the magnetosphere 

We need to know a soft photon field to calculate the pair-creation rate rj p ±. For GeV 7-ray 
photons, which dominate in the curvature photons radiated in the gap, a p is maximized at 
the soft X-ray band. The observed X-rays are dominated by the surface blackbody radiation 
for the Vela-like pulsars, and by the non-thermal radiation for the Crab-like pulsars. In this 
paper, because we apply our model to the Vela-like pulsars (JQ), we take only the surface 
blackbody radiation as the soft photon field. 

For the blackbody case, the photon number density between energy m e c 2 ex and m e c 2 (ex + 
dex) is given by the Planck law, 



where r is the radial distance from the centre of the star, A s is the observed emitting area 
and kT s refers to the surface temperature. We calculate the pair-creation rate (JTHJ) by using 
(1231) at each point in the gap, where the collision angle /i c is given by 



3 CURVATURE RADIATION FROM THE OUTSIDE OF THE GAP 

In the HS model, a self-consistent determination of the accelerating electric field is successful 
in reproducing the 7-ray spectrum in the GeV band, at which the pulsed radiation peaks 
in vF v diagram. However, the calculated spectrum deviates from the observed that in lower 
energy bands. This problem is solved by treating properly how particle energy evolve near 
and beyond the gap boundaries. 

Because the saturated Lorentz factor r sat given by (jlfij) is proportional to E^ 4 , T sat 
does not change much in the gap, except for around boundaries. In such case, the curvature 
spectrum produced in the gap is similar to the spectrum that is emitted by mono-energetic 




(23) 



jj, c = cos0 7 x sin#, 



(24) 



with 7 x = sin x (r sin 9Q/c). 
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particles, and follows E~ p with the photon index p ~ 2/3 below the cut-off energy E c . 

However, observations indicate much soft slopes; E^ 1 - 3 — E~ 1S in MeV band. 

If one uses the saturated Lorentz factor given by equation (|16|). T sat decreases near the 

inner and outer boundaries because E\\ decreases. It is notable that T sat is the Lorentz factor 

at which t ac = td- The time scales t ac and become longer near the boundaries as compared 

with the crossing time t cr . In fact, the ratio of t cr to t d , 
t d 3 m e c 2 R 2 C . 

~ g 3 x 

ai^e / n y 1 (i^y 3 / R c \ 2 

Wjl UOOrads- 1 / UO 7 / V0.5^J ' 1 ' 

is larger than unity near the boundaries. As a result, the Lorentz factor is kept greater 

than r sa t near and beyond the boundaries. These particles escape form the gap with a large 

Lorentz factor T out ~ a few times 10 7 and will radiate 7-rays even outside of the gap. The 

damping length after the escape should be comparable with or longer than the gap width 

W\\. 

On these aspects, we take account of the radiation from the outside of the gap as well 

as that from inside. Although Hirotani et al. (2003) calculated T out numerically by solving 

the unsaturated motion of particles in the gap, we shall estimate T out semi-analytically; 

we assume that the particles are ejected from the gap with the Lorentz factor given by 

a saturated value at t cr = t d . On leaving the gap, particles simply lose their energy by 

radiation, because E\\ is screened out. The Lorentz factor of the particles ejected from the 

gap decreases as 
dT(s) = 2 e 2 r 4 ( S ) 

ds 3m e c 2 R 2 c (s)' { ' 

The rate of decrease in V per unit length is proportional to T 4 (s) if i? c =const. Then, the 

spectrum of 7-rays from the outside of the gap is the same as the spectrum made by particles 

with an energy distribution of 

^MQ K r-*. (27) 

In short, we obtain a power- low 7-ray spectrum by superposing the curvature radiation at 
different positions; therefore, we need not assume a power energy distribution for particles 
to explain the observed power spectrum. 

The curvature photons emitted by the particles obeying equation (|27jl exhibits a soft, 
power- low spectrum with a photon index around p = 5/3 for E^ < E out where E out = 
(3h/4n)(Tl ut c/ R c ) is the critical energy for the highest-energy particles. The soft, power-low 
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spectrum will extend down to the energy E min = (3/i/47r)(r^ n c/i? c ), where T min is the 
Lorentz factor below which the 7-rays are beamed out from the line of the sight due to the 
magnetic curvature, for instance. Below E m i n , the curvature spectrum becomes hard with 
a photon index around p — 2/3. If we estimate E min with the particles energy when they 
escape across the light cylinder, we obtain E min ~ 80 MeV for the Vela pulsar. For example, 
the outwardly moving particles escape from the light cylinder with T ~ 8 • 10 6 (for the Vela 
parameters). Thus, we expect that the calculated spectrum is improved between E min and 

Eout- 

The 7-rays radiated towards the star may be converted into e^ 1 pairs by 7B — > e ± B 
process and cause the pair-creation cascade. Cheng & Zhang (1999) showed that the syn- 
chrotron spectrum from these cascade pairs has the photon index of p ~ 1.9 in X-ray region. 
This spectral index is consistent with the pulsed RXTE soft component for the Vela pulsar 
(Harding et al. 2002). According to the model of Cheng et al. (2000), the 7-rays escaping 
from the gap through the inner boundary are much fainter than 7-rays through the outer 
boundary because the pair-creation occurs more efficiently near the inner boundary. In the 
present paper, the luminosity of the inwardly and the outwardly propagating 7-rays depend 
on the assumed currents j\ and j'2, respectively. In §H we will find ji ^> j'2 to reproduce the 
observations. With this parameters, we have much more outwardly propagating 7-rays than 
inwardly. 

The curvature of field lines will be different from that of dipole near the light cylinder. 
In this paper, therefore, we investigate the effect of the non-dipole field on the curvature 
spectrum; this effect was not considered in Hirotani et al. (2003). The many authors (e.g., 
Michel 1973; Contopoulos, Kazanas Sz Fendt 1999) have attempted to set the structure of 
the field lines, but it is yet far from an enough understanding of the actual structure of the 
magnetic field (e.g., non-aligned rotator with the pulsar wind). To examine the variation 
of the curvature radius of field line near to the light cylinder, we calculate each curvature 
spectrum for the following three configurations: 

Case 1: the field line is defined by a pure dipole field that is tangent to the light cylinder 
(so-called the last-open field line for the dipole field); 

Case 2: the field line is the same as Case 1 below w = 0.6wi c , but deviates beyond 0.6zui c 
so that the radius of curvature may change as R c oc 1/(1 — uj/wi c ); 

Case 3: Below 0.Qzui c , the field line is defined by a pure dipole field that is tangent to the 
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4.4 0.5 0.6 0.7 0.8 0.9 1 

distance from the rotation axis (light radius) 

Figure 1. The curvature radius of last-open field line for three cases when oj;„=45°. The solid, dashed-line and dotted-linc 
correspond to the Case 1 (the case for pure dipole field), Case 2, and Case 3, respectively. 



cylinder of radius 0.8wi c . Beyond O.Qwi c , the field line deviates from the pure dipole, and 
the radius of curvature of the field changes as R c oc 1/(1 — zu/-aj lc ). 

Case 2 and Case 3 simulate the effect of the increase of open magnetic fluxes. In Case 3, 
we take account of a possibility that the last-open field line changes inward dipole field line 
as compared with that of Case 1 by the influence of the pulsar wind or the electric current; 
for instance, the colatitude of the foot-point for the non-inclined pure dipole field line on 
the stellar surface is Of ~ (R^/zuic) 1 ^ 2 for Case 1 (for Case 2) and Of ~ (R^/O.Szuic) 1 ^ 2 for 
Case 3, where R* is the star's radius. For the axisymmetric force-free magnetosphere, it is 
indicated that the field lines around the last-open field line deviate from dipole to radial 
configuration from w ~ 0.6tu; c , and have a point where R c = oo near the light cylinder 
(e.g., fig. 3 in Contopoulos et al. 1999). These effects are simulated in Case 2 and Case 3. 
Fig- CD shows the curvature field for Case 1 (solid-line), Case 2 (dashed-line), and Case 3 
(dotted-line) between w = 0Azui c and the light cylinder when the magnetic inclination cti n 
is 45°. 

The electron-positron pairs produced by pair-creation will radiate by synchrotron process 
in the outside of the gap. These pairs are produced with significant pitch angle Xi an d 
T ~ 5 • 10 3 . The critical photon energy of the synchrotron radiation around the outer-gap 
becomes 

3ehT 2 B 
E c = sin x 

A7rm e c 



4.4- 10 5 



(f^H^gH* (28) 



It is, therefore, entirely fair to neglect the synchrotron component emitted the outside 
of the gap when we consider a 7-ray spectrum above 100 MeV. 
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Table 1. The observed parameters 



Pulsar 


d 


n 


Bl2 




A s 




kpc 


rad/s 


10 12 G 


eV 


km 2 


Vela 


0.5 


70.6 


3.4 


150 


6.6(d/0.5kpc) 2 


B 1706-44 


1.8 (DM)/2.5 (H I) 


61.6 


3.1 


143 


13(d/2.5kpc) 2 



d is the distance from the earth to the pulsar. For PSR B1706-44, we adopt the 1.8 kpc 
proposed by dispersion measure and 2.5 kpc by H I absorption as the d. B12 represents 
the magnetic field strength at the star surface in units of 10 12 G. The X-ray data are 
referred from Ogelman, Finley & Zimmerman (1993) for the Vela pulsar and Gotthelf et 
al. (2002) for PSR B1706-44. 

4 APPLICATION 

In this section we apply our model to the twin pulsars, PSR B0833-45 (Vela) and Bl 706-44, 
whose observed parameters are listed in Table [TJ The model parameters are (jtat, ji, 32), 
the magnetic inclination angle «j n between the axes of rotation and magnetization, and 
the cross section area A gap of the gap. The model parameters (jtot, ji, 32) are chosen so as 
to reproduce the observed spectral properties. We adjust A gap , which should not exceed 
(zui c ) 2 , to the observed flux. For PSR B1706-44, the proposed distance to the pulsar is about 
1.8 kpc by the dispersion measure and about 2.5 kpc by the H I absorption. The X-rays 
emission from PSR B1706-44 were reported for the first time by Gotthelf et al. (2002) with 
the Chandra X-ray observatory. 

4.1 The electric structure 

In Fig. |2l we present the calculated En for Case 1 for both pulsars. The model parameters 
are {jtati ji, 32) — (0.201,0.191,0.001) and a in = 45°. The abscissa refers to the arc length 
from the pulsar surface along the last-open field line. 

We see that the gap width Wu is shorter than ci7j c , which is 4.25 x 10 8 cm for the Vela 
pulsar and 4.86 x 10 8 cm for PSR B1706-44. In the present model, W\\ is characterized by 
the mean free path for the pair-creation as Appendix [X] shows. Since the gap is filled with 
abundant 7-rays and X-rays, the mean-free path is shorter than w\ c . The ratio of the mean- 
free path for the Vela pulsar and that for PSR B1706-44 becomes ~ 0.9. As a result, we 
obtain the shorter gap width for the Vela pulsar (solid-line) than that for PSR B1706-44. 

As Fig. 121 shows, a greater En is obtained for the Vela pulsar than for PSR B1706-44. 
This can be understood from the fact that the strength of En increase with W\\ and with the 
magnitude of pcj in the gap. The ratio of pcj for the Vela pulsar and that for PSR B1706- 
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8000 



Vela — 
B1 706-44 (1.8 kpc) ■ ■ ■ 
(2.5 kpc) — 



q 



1.6 



1.7 



1.8 1.9 2 2.1 2.2 2.3 

distance along field line (10 A 8 cm) 



2.4 



2.5 



Figure 2. The acceleration field for the Vela pulsar (solid-line) and PSR B1706-44 (dotted-line for d = 1.8 kpc and dashed-line 
for d = 2.5 kpc). (jtot, 31, J2)=(0.201, 0.191, 0.001) and a in = 45° 

44 is ~ 2. The effect of the difference in pcj overcomes that in the gap width, that is, 
PSR B 1706-44 has a greater gap width W\\ than the Vela pulsar does. 

Fig. Elalso shows how E\\ depends on the assumed distance for B1706-44. The gap width 
for d = 1.8 kpc (dotted-line) is wider than that for d = 2.5 kpc (dashed-line) because the 
X-ray luminosity decreases with decreasing assumed distance. The decreased X-ray photon 
density results in an increased pair-production mean-free path, and hence increased gap 
width. Therefore, if we adopt a nearer distance to the pulsar, the calculated En becomes 
large as Fig. |21 shows. This also causes a harder spectrum in the case of d = 1.8 kpc than 
the case of d = 2.5 kpc (Fig. @J. 

4.2 The 7-ray spectrum 

The dimensionless current parameters are chosen so as to reproduce the observation. We 
obtain the best-fitting parameters (jtot, ji, jz) — (0.201,0.191,0.001) for the Vela pulsar 
with ai n = 45°. Fig. 01 shows the best-fitted curvature spectrum (solid-line) for outwardly 
propagating 7-rays, which is composed of radiation from the inside (dashed-line) and the 
outside (dotted-line) of the gap for Case 1. For comparison, we plot the phase averaged 
EGRET spectrum. It's peak and cut-off energies are in good agreement with the EGRET 
observations. 

As we expected in the spectrum from the outside of the gap appears between E out ~ 
1 GeV and E min ~ 80 MeV with the photon index p ~ 5/3, where E out and E min correspond 
to the critical energy of curvature radiation of the particle's Lorentz factor at the boundary 
and the light cylinder, respectively. Below E min , the spectrum has p ~ 2/3. By including 
the emission from outside of the gap, we can improve the spectrum above 100 MeV. 
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(1994) for EGRET and Strickman et al. (1996) for OSSE and COMPEL. 



It follows from Fig.|3]that the calculated total spectrum in 100 MeV-GeV bands becomes 
slightly harder than the EGRET data. This is because we assume that the particle's motion 
is saturated at the equilibrium Lorentz factor ()16|). which overestimates the particle energy, 
and hence the 7-ray luminosity. A better predicted spectrum could be obtained if we would 
consider a unsaturated motion for the particles as Hirotani et al. (2003) did. Nevertheless 
the outside gap emission can also be treated properly with the Lorentz factor solved with 
equation (j2T?j) . 

Fig. E] shows the total spectrum for PSR B1706-44. It is noteworthy that the same 
dimensionless parameters (jtot, 31,32), and Oj n as the Vela pulsar reproduce PSR B1706-44 
data as well. However the observed peak energies for the two pulsars are different: About 3 
GeV for the Vela pulsar, and about 1 GeV for PSR B1706-44. The difference is mainly due 
to the rotation periods. As showed in ^4.11 the difference in pcj for the two pulsars appears 
as the difference of -Ej|. As a result, the saturated Lorentz factor given by equation (|16|) and 
the obtained spectral peak energy for PSR B1706-44 are less than those for the Vela pulsar. 
Thus, we need not change dimensionless parameters (jtot, ji, 37) • 

We fit the flux, by changing the cross section of the gap, A gap . We use A gap = (5W\\) 2 ~ 
(0.5w lc ) 2 for the Vela pulsar, and A gap = (10VK M ) 2 ~ wf c for PSR B1706-44. The value 
of Ag a p for B1706-44 appears to be too large, since the gap is located at s ~ 0.5zui c . In 
fact, we have A 9ap /47r(0.5^ c ) 2 ~ 1/3 for B1706-44, which means that the hemisphere are 
almost covered with the gap. Such a large A gap may be unrealistic. This flux problem will 
be discussed in £0 

Let us examine how the spectrum depends on the magnetic inclination, a« n . Romani & 
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Yadigaroglu (1995) estimated the inclination angle to be a in ~ 65° for the Vela pulsar and 
ai n ~ 45° for PSR B1706-44 to explain the phase difference between the radio and 7-ray 
pulses. We also show the calculated spectrum in the case ai n = 65° for the Vela pulsar in 
Fig. El where we adopt (j tot , j x , j%) = (0.202,0.192,0.002) to explain the observed spectral 
peak energy. The strength of magnetic field in the gap increase with oii n because the crossing 
point of the null surface and the last-open field line comes close to star. On the other hand, 
the gap width is shortened because the increased soft photon density results in the decreased 
pair-creation mean-free path for the inwardly propagating 7-rays. The former effect tends to 
increase E\\, while the latter to decrease it; therefore, the two effects compensate each other. 
As a result, we obtain the observed spectral peak energy by using almost the same current 
parameters for = 45°. On these ground, we suggest that almost the same currents in 
unit of the GJ value are running through the gap for the Vela pulsar and PSR Bl 706-44. 
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4.3 Spectral dependence on the curvature of field line 

Let us consider other configurations of the field line. Fig. El represents how the particle's 
Lorentz factor changes along the field lines with different field line configurations, where the 
Vela parameters are used. Three cases, Case 1 (solid), Case 2 (dashed) and Case 3 (dotted), 
are examined as had been described in JJJ The Lorentz factors for Case 2 and Case 3 change 
little between w = 0.6zui c and w\ c because of the stretched magnetic field lines. 

Fig. [7| shows the calculated total spectra for the three cases. The dependence of each 
configuration of the field line appears in the spectra around 100 MeV. On such a stretched 
field line as in Case 2 or Case 3, curvature process little contributes to the spectrum. Because 
of the less efficient curvature emission, particles escape across the light cylinder with a large 
energy as compared with a pure dipole field. As a result, we obtain a greater lower-cut off 
energy E min ~ 200 MeV for Case 2 (cf. E min ~ 80 MeV for Case 1). 

This effect is expected to appear in the observed spectra. Therefore, it is possible to di- 
agnose at some level the actual field configuration in the magnetosphere from the observed 
spectra. For example, the EGRET, COMPTEL and OSSE observations indicate that the 
spectrum turns over below 100 MeV, that is, E min < 100 MeV. Possibly a dipole-like field 
is preferable near the light cylinder than a stretched field, which is expected from magneto- 
hydrodynamics. 

The effect of stretched fields will be manifest in a phase-resolved spectrum rather than a 
phase-averaged one. According to the double-peaked 7-ray models (e.g., Romani 1996; Cheng 
et al. 2000), the 7-ray photons radiated near the light cylinder contribute to the first pulse- 
phase, which leads the bridge-phase. Therefore, the stretching effect is expected to appear 
in the spectrum of the first pulse-phase. A diagnosis of the field line curvature is, therefore, 
carried out simultaneously with a position-pulse-phase mapping, which is necessary when we 
model the gap in a three-dimensional space. Thus, one can examine the field line geometry 
by using the fact that the curvature spectrum around 100 MeV depends on the field line 
geometry. On the other hand, in the phase-averaged spectrum, the radiation on stretching 
fields may be buried in the curvature radiation of low energy particles near the upper surface 
in the gap, where En is small owing to the screening effect. 

Even if a possible variation of the field line curvature is taken into account, the discrep- 
ancy in the COMPTEL and OSSE bands cannot be resolved. It seems difficult to diagnose 
the field line with these observations. In general, if the curvature radius varies smoothly on 
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the particle's path to the light radius, we are able to think that the particles escape from 
the light cylinder with the Lorentz factor at which ctd ~ toi c , 
fi X- 1 / 3 f<R c >\ 2 ^ 



r *~ 117 x 10 Us) (-^-) 



(29) 



where < R c > is the averaged curvature radius of the field line along which the particles 
move. For this Lorentz factor and the radius of curvature, particles emit 7-rays at energies 
around 



E l : ^ 160 ( K Rc> ) MeV. 



(30) 



which may appear in the curvature spectrum as the lower energy cut-off E m i n . If the ob- 
served break at ~ 10 MeV for the Vela pulsar were caused by this escape effect, we would 
have < R c >~ 0.1wi c , which should be too short. On these ground, we conclude that the 
curvature radiation from the primary particles does not account for the COMPTEL and 
OSSE components. 
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5 DISCUSSION 

In short, we showed that curvature flux emitted outside of the gap dominates that emitted 
within the gap around 100 MeV. With this correction, the model spectrum was improved 
above 100 MeV. We then found that the field line curvature near the light cylinder affects the 
curvature spectrum in several hundreds MeV. We also showed that the difference in spectral 
peak energies of the Vela pulsar and PSR B1706-44 is explained by intrinsic difference in the 
magnetospheric parameters of the individual pulsars. We need not change the dimensionless 
current parameters. 

5.1 Dominance of ji to total current 

In £0J we find that j 2 should be much smaller than jx- In this subsection, we consider the case 
in which j 2 is not smaller than j\. For a large inclination, (e.g., for o%j greater than or nearly 
equal to 45°), the gap width W\\ is mainly determined by the pair-creation mean-free path 
for the inwardly propagating 7-rays because the collision angle with X-rays is nearly head- 
on for the inwardly propagating 7-rays and nearly tail-on for the outwardly propagating 
7-rays. In this case, W\\ sharply decreases with increase of j 2 . Because E\\ decreases with 
Wu, the spectrum becomes soft by increasing j 2 . Fig. |H] shows the spectra for j\ = 0.191 and 
j g = 0.009, and j 2 = 0.001 indicated by solid-line, but j 2 = 0.021 by dashed-line for the Vela 
pulsar. A choice of j 2 = 0.021, which is significantly larger than the best-fitting value but 
smaller than j\, does not reproduce the observed spectrum. 

Although the cases of ji < 0.1 can reproduce the observed spectral peak by appropriately 
choosing j 2 , it is difficult to reproduce the observed flux with a cross section that is smaller 
than (zuic) 2 . Therefore, a large ji and a small j 2 are preferable. 

5.2 Comparison with previous models 

Romani (1996) discussed the 7-ray spectrum for the Vela pulsar by curvature radiation, 
assuming a specific energy distribution dN e /dE\\ oc Eu to connect the assumed electric 
field, En oc 1/r, and the number density of the particles at the upper surface of the gap. 
This relation gives the particle energy spectrum dN e /dT oc T™ 4 , which is identical with 
equation (|2*Tj) . as long as the particles have the saturated Lorentz factor ()16j) . The curvature 
radiation of the particles produced in the gap dominates in the spectrum (see, fig. 5 in Romani 
1996). 
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On the other hand, in the present model we obtained equation (|27j) as a result of the 
radiative loss of the particles in the outside of the gap. The pair flux produced in the gap 
is little, j g = 0.009 for instance; therefore, the injected particles j\ play a crucial role in the 
formation of the curvature spectrum above 100 MeV. 

This difference in the origin of the 7-rays between the two models will come from the 
difference in geometrical picture of the gap, slab-like outer-gap model (Romani's picture) or 
thick gap model (present model). The slab-like gap model, the pioneering studies for which 
were done by Cheng et al. (1986a,b), assumes that the gap width along the magnetic field 
line is much larger than the trans-field thickness of the gap on the meridional plane. In 
this slab geometry, the number density of pairs near the upper surface of the gap grows 
exponentially in the trans-field direction by the pair-creation, because the 7-rays radiated in 
the gap can effectively across the field lines before materializing as pairs. In the pair-creation 
region, because E\\ is not abruptly screened, a copious created particles are accelerated and 
emit the observable 7-rays. 

In the present model, on the other hand, we assume that the gap width is much shorter 
than trans-field thickness. In this case, the radiated 7-rays do not effectively across the 
field lines in the gap before materializing; therefore, we can consider the gap with the one- 
dimensional model. The calculated E\\ in this non- vacuum model decreases with increasing 
of pairs produced in the gap because of the screening effect. As a result, the cases of the large 
pair-flux production in the gap, the large j gap in (I22|) . do not reproduce the observed spectral 
peak around GeV. As discussed in ^5.11 therefore, we give a large injected particle flux into 
the gap through a boundary with a small j gap to explain the observations. If we consider 
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the Romani's picture within our electrodynamical treatment, we should deal with a model 
of two- or three-dimensional geometry because the present one-dimensional assumption can 
not be used for such slab-like geometry. 

5.3 Three dimensional electrodynamical model 

Although the difference in spectral shape of the twin pulsars, the Vela pulsar and PSR B1706- 
44, is explained with the same dimensionless current parameters, the difference in the flux 
is difficult to be understood. Because the two pulsars have almost the same magnetospheric 
parameters and the surface temperature, one may expect similar luminosities. If intrinsic 
luminosity were the same, difference in flux would be attribute to the distance so that the 
expected ratio of flux may be between 13 and 25. However, the observed value is about 5. 
Zhang and Cheng (1997) estimated the transverse thickness of the outer-gap by considering 
a condition for 7-rays to materialize marginally in photon-photon collisions, which may 
explain difference in A gap between the two pulsars. They obtained 

/ ~ 5.5P 26 / 21 B u i/T (31) 

where / is the transverse thickness in unit of the light radius, P is the rotational period. 
With this weak dependence on the magnetospheric parameters, one cannot expect difference 
in A gap between the two pulsars. One possibility is difference in viewing angle. For this point, 
we cannot proceeds without extending the present model to a three-dimensional one. We 
need to solve the trans-field structure of the electric field and the current distribution. 

There is another limitation of the present one-dimensional approach. To get a reasonable 
spectrum, the dimensionless current is should be less than or nearly equal to 0.2 in units of 
the GJ value. Within the frame work of the one-dimensional model, this value cannot be 
enlarged because the gap is quenched by large currents. This fact limits the total luminosity 
seriously. On the other hand, a current of the order of the GJ value should circulate in the 
pulsar wind, and it is the current that runs through the gap on its way to the star and 
contributes to the spin-down. In the outer-gap model by Cheng et al. (1986a), pairs are 
created much in the convex side of the gap, and a large current is flowing there. To evaluate 
the net current in the gap and the resultant luminosity correctly, the trans-field structure 
must be taken into account. 

Finally, we note that the present model can only be compared with the observed phase- 
averaged spectra. On the other hand, the phase-resolved spectra have been discussed with 
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three dimensional geometry such as Romani (1996) for the Vela pulsar and Cheng et al. 
(2000) for the Crab pulsar. We could discuss this issue from the electro dynamical point of 
view if we would extend the present approach into a three-dimensional model 
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APPENDIX A: THE ELECTRIC FIELD SCREENING AND THE GAP 
WIDTH 

We give the gap width W\\ by solving the continuity equations (JOJ), (fT2"|) with some approxi- 
mations. 

We integrate the Poisson equation (J3j) along the field line from s x to S2 to obtain 

S2 [p(s)-p GJ (s)}ds = 0, (Al) 
1 

namely, the total amount of the effective charge in the gap equals zero. For the case of the 

vacuum gap (p = 0), the null point, pcj = 0, is located between s% and 82- For a given si, 

one can find S2 so that (|A1J) may be satisfied. The gap width is arbitrary. 

For a non- vacuum gap, Wu is characterized by pair-creation mean-free path. Since the 

gap width is shorter than w\ c (JQ), we assume the physical quantities in the gap to be 

constant value except for the particle's number densities (N±) and the distribution function 

(G±) for 7-rays. Then equations © and (fT2"J) are combined to give 
d 2 N±(s) 



ds 2 
where 



±[a + N + (s) - a_JV_(s)], (A2) 



a± = / de^r] p± (e^)r] c (e^). (A3) 
cr Jo 

With this solutions for N±(s) and the boundary conditions (J2DJ), (EH), we obtain the following 
representation for the gap: 



JO 



1/0;+ + a_ sinh( A /a + + oZWn) 

I + \ JL a++a-J \ JZ a++a-J _ g 

tanh( v /a + + aIW\\) ' 
For a largely inclined pulsar (e.g. for oti n greater than or nearly equal to 45°), inwardly 
propagating 7-rays collide almost head-on with the surface blackbody X-rays. In contrast, 
outwardly propagating 7-rays collide almost tail-on. Due to this difference, we obtain 

w ''~(7^ cosh i^)- (A5 » 

This equation indicates that W\\ for a large inclination is characterized by the pair-creation 
mean-free path of inwardly propagating 7-rays unless j'2 = 0. 

By giving (jtoti hi 32), the position of the gap in the pulsar magnetosphere is determined 
by the set of S\ and S2 that satisfies relation (jA4j) and condition (jAljl . Thus, it does not 
need to have the null point in the gap. 



